Effect of weak disorder induced on the first order transition in a p-wave 

superconductor. 
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We investigate the influence of weak disorder on the fluctuation driven first order transition 
found in a p-wave superconductors. The flow equations for the disorder averaged effective theory is 
derived in an 4 — e expansion. The fixed point structure of these flow equations and their stability 
is discussed. 



The role of quenched disorder near second order phase 
transitions have been the subject of intense investigations 
in the last few decades. These investigations led to the 
development of the so-called Harris criterion [B| which 
helps to gauge the relevancy of weak disorder near the 
clean critical point. Also, there is a significant body of 
work that deals with strong disorder effects near various 
critical points. These are generically termed as Griffiths 
phenomenon which leads to very interesting behav- 
ior of various observable quantities like the susceptibility 
in the vicinity of critical points. However, the role played 
by quenched impurities near first order transitions has 
not been so thoroughly investigated. A notable excep- 
tion to the above statement is the pioneering work of 
Aizenmann and Wehr, [Tfl where it was proved that in di- 
mensions, d < 2 the addition of quenched disorder leads 
to rounding of the transition. Hui and Berker, ||, also 
came up with exactly the same conclusion using heuristic 
Renormalization Group (RG) arguments. 

The case in higher dimensions is not so clear. In other 
words, in certain scenarios, the presence of quenched dis- 
order rounds the first order transition and in other cases 
it seems that there is no impact on the order of the tran- 
sition. For instance, the first order transition inherent in 
a magnetic system with cubic anisotropy can be shown 
to be unaffected by the addition of weak-disorder [Io| . 

An example in higher dimension wherein disorder does 
alter the nature of the transition is given by the case of 
type-II superconductors Q : In the case of a clean type-II 
superconductor Halperin and co-workers had shown that 
the transition was of the fluctuation driven first order 
type. This conclusion was based on a 4 — e expansion, 
wherein the RG flows showed runaway behavior. These 
results for the clean superconductor were also in accor- 
dance with the results obtained by Coleman and Wien- 
berg [jj for the same model. In Ref. 0] , Boyanovsky and 
Cardy analyzed the effect of impurities on the fluctuation 
driven first order transition in such a s-wave supercon- 
ductor. In particular it was shown in Ref. [2], that if one 
adds quenched disorder as a random mass term, then as 
long as the "bare disorder" strength is above a threshold 
value the first order transition inherent in scalar electro- 
dynamics, (type-II superconductor), gets converted into 
a second order transition. In other words, if one starts 



with a bare disorder value that is below the threshold 
value then the first order transition persists. These re- 
sults were obtained in an e = 4 — d expansion. 

In similar vein, using the e = 4 — d expansion, Bla- 
goeva, et al., [HI], studied the nature of the transition 
in an unconventional superconductor with orthorhom- 
bic and cubic anisotropics: They found that the phase 
transition in these superconductors was predominantly 
weakly first order. However, in a situation where the or- 
thorhombic anisotropy is absent, they found a new crit- 
ical fixed point. In [17| . the influence of weak random- 
ness on the fluctuation driven first order transition in 
such unconventional superconductors were studied. They 
concluded that in an e = 4 — d expansion the transition 
in such anisotropic unconventional superconductors re- 
mained weakly first order. 

Thus, from the discussions of the preceding paragraphs 
it can be seen that in d > 2, the influence of quenched 
disorder on first order transitions is quite model depen- 
dent. 

In this brief communication, we look at the effect of 
quenched non-magnetic impurities on the critical behav- 
ior of a p-wave superconductor: The clean (impurity- 
free), version of this model was recently studied by Qi 
Li and collaborators 0]. In an e = 4 — d expansion, 
when the number of order parameter components, to, are 
small, they showed that the phase transition in such an 
unconventional superconductor is also of the fluctuation 
driven first order variety. These results were analogous 
to the results obtained for the case of the s-wave super- 
conductor that we had elucidated earlier. To in- 
corporate the effect of disorder, we start with a random 
mass generalization of the effective action derived by Qi 
Li and collaborators for the p-wave superconductor. Un- 
der this generalization, we set the distance to criticality 
t—*t — S(t(x)), and we obtain: 
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Here, D = V — i2eA, is the minimal coupling term. A is 
the vector potential, and if>, is the m component complex 
order parameter field. We integrate out the quenched 
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disorder by using the replica trick. The replica trick 
essentially results in re- writing the F = — {\ogZ} d{s = 



lim„^o 



{^"}d 



Here, the subscript just implies that 



one has to perform a disorder average. Now, for cal- 
culational ease we assume that the disorder distribution 
is delta correlated. In other words, we assume that the 
{8t(x)}dis = 0, and {St(x)6t(y)} < a s = A6(x — y). We now 
perform the disorder average and then perform a 1-loop 
RG calculation. The result of which this analysis gives 
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The flow equations, Eq. [2j appropriately rescaled, 
maps onto the case of the disordered s-wave supercon- 
ductor studied in Ref. 0] when v = 0. In the same 
spirit, it is identical to the flow equations derived in Ref. 

for the p-wave superconductor when A = 0. In what 
follows, we analyze the fixed points of the above set of 
equations. We will linearize the flows around the fixed 
points and thus discuss their stability. Throughout the 
course of this paper we have suppressed the flow of the 
distance to criticality t. This is on account of the fact 
that the t is always a relevant operator. Thus, we look at 
only the flow in the (u,v,A) space. Furthermore, since 
we are interested in the superconducting transition, we 
will look for a fixed point with finite value of e 2 . 
Case 1 : u* = v* — A* = 0. This is the Gaussian fixed 
point which has eigenvalues A„ = A„ = Aa = 4 — d. This 
implies that the Gaussian fixed point is unstable for all 
d < 4, consistent with the naive power-counting. 

Next we look at a class of fixed points where one of 
the coupling constants is constrained to a non-zero fixed 
point value. Now, the three such fixed points are : 

Case 2 : u* ^ 0, v* = 0, and A* = 0. 
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Consistent with Refs. 0], 0], we call this pair of fixed 
points the pure s-wave fixed points. It can be shown 
that this fixed point is physical, (i.e., u* is real), only 
when m > m c w 183. It can be shown that for physi- 
cally relevant parameter values, i.e., m > m c , this pair 
of fixed points are unstable. The most obvious route 



to show the instability of the pure s-wave fixed point is 
to calculate the eigenvalue A„, which can be shown to 
be: A„ = ^r(2(rn - 2)u* 2 + ^pr). As can be seen this 
eigenvalue is always positive in the physically relevant 
parameter regime. 
Case 3: u* = 0, v* ^ and A* = 



u* = 
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A* = 



Once, again this is an unstable fixed point. It can be 
shown that one of the eigenvalues, namely Aa = — 2(m — 
4)u*, is positive for all m > 4. Now, for m < 4, even 
though Aa switches signs, it can be shown by following 
Ref. [3J], that this fixed point is unstable. 
Case 4: u* = 0, v* = and A* ^ 
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This is an un-physical fixed point, as A*, which is related 
to the width of the disorder distribution, is negative at 
this fixed point. Since, this is an un-physical fixed point, 
we will refrain from discussing its stability here. 

Now, we turn our attention to the case where two of 
the three coupling constants attains non-zero fixed point 
values: 

Case 5: u* ^ 0, v* = and A* ^ 
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Since the fixed point value v* = 0, we will call this case 
the disordered s-wave fixed point. It is seen from follow- 
ing Ref. 0, that the A u and Aa are both negative for 
small number of order parameter components. Thus, we 
need to study the eigenvalue A„ that we obtain from lin- 
earizing around the v* — fixed point. It can be shown 
that At, = (4(m — 2)m*+4A*). It is obvious that only for 
the case m = 1, will this eigenvalue be negative. Thus, 
implying that this fixed point is a critical fixed point for 
m = 1. For all other values of m, this fixed point is un- 
stable and the flow is away from this fixed point. It is of 
little surprise that at the disordered s-wave fixed point, 
the operator v is irrelevant, as it can be shown that the 
v term in the action is identically zero. Thus, for m = 1, 
we flow into a critical fixed point first obtained in Ref. 
0. As can be seen from Fig. [TJ there is a minimum 
critical value of the bare disorder strength beyond which 
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Case 6: u* ^ 0, v* ^ and A* = 



u -A plot with m=1, v=0.05 and e=0.1 
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FIG. 1: Numerical solutions of the RG equations in the u — A 
plane. The plot is for the case m = 1, the only case where 
there is a fixed point for low enough number of order param- 
eter components. It is clearly shown that when there is a 
sepratrix in this set of equations. In other words, when the 
bare disorder is less than some threshold value then one flows 
away from the fixed point. However, when the disorder is 
greater than the threshold value, we flow into the fixed point. 
These are consistent with the results of Q]. 




6 



c 

12 



m 



v 

A* 



e [(to - 3) (to + 4) (to + 18) ± 3Mi 



2 TO 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Where, Mi = (m 4 - 204to 3 - 1356to 2 - 864m - 15552) l / 2 
and D\ = to 3 + 4m 2 — 24m + 144. It can be shown by 
following the arguments in that such types of fixed 
points are only physical when to > m' c . Here, m! c ~ 
210.45. 

Furthermore as argued in Q we can show that when 
to > m' c , we get a pair of such fixed points with the 
one with the higher numerical value of v* and u* , being 
the one with negative eigenvalues \ u and At,. Since, the 
fixed point with the lower numerical value of u* and v* , 
have been shown to be manifestly unstable in Q , we will 
ignore it from our future discussions. In analogy with 0], 
the fixed point with larger value of u and v will be called 
the pure p-wave fixed point. 

Now, all that remains to be seen is whether this pure 
p-wave fixed point is stable to infinitesimal perturbation 
in A. This can be easily accomplished by evaluating 
Aa, which has the form 4(2 - m)u* + 2(4 - m)v* . By 
substituting the functional form for the fixed point into 
the expression for Aa , we can easily show that that Aa is 
negative for to > m' c . These expressions are very cumber- 
some and thus we will discuss our results for particular 
number of order parameter components. For instance, 
for the case to = 300, we can show that the eigenvalues 
are: -1.04576e, -0.537164e and -0.52356e which clearly 
indicates that the fixed point is stable. Similarily when 
to — » oo, it can be shown that the eigenvalues are all 
negative at the clean p-wave fixed point. Thus, the pure 
p-wave fixed point is stable for large number of order 
parameter components. 

Case 7: u* = 0, v* ^ and A* ^ 



FIG. 2: Flow plot from RG equations for two different values 
of m, m = 1 and m = 2 . Both flows start from the same 
initial value of u = 0.5, v = 0.04 and A = 0.1. For m — 1, 
the flow is into the fixed point given by u* — 0.675 and A* = 
0.4375. But for m = 2, as the graph clearly shows there is 
runaway flow towards large values of u and A. In fact as seen 
in the inset for m = 2, the flow crosses over into the classically 
unstable regime, u < 0. 



the system goes through a second-order transition. Else, 
the transition remains first order. 
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As can be seen from the above set of equations, the fixed 
point value A* is positive only for m = 5. Thus, it is 
sufficient to discuss the stability at this particular value 
of m. For m = 5, it can be shown that the eigenvalues are 
29.6994e, -7.94968e± 4.77846d and -e. Thus, implying 
that the fixed point is once again unstable. 
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Case 8: u* ^ 0, v* ^ and A* ^ 
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where M 2 = (to 4 + 1740to 3 - 48012to 2 
15552) 1 / 2 and D 2 = (2m 3 - 55m 2 + 114to - 36) 

These pair of fixed points will be called the disordered 
p-wave fixed points. As can be seen the expressions for 
the fixed point values are extremely cumbersome. So 
is also the case for the eigenvalues that we obtain as 
a result of linearizing around this fixed point. Thus, 
for the sake of simplicity, we will discuss the behav- 
ior of the system near these fixed points for particular 
representative values for the number of order parame- 
ter components. For small values of m, we can show 
that the fixed point is unphysical. For example, for 
to = 3, it can be shown that the fixed point value is u* — 
0.291667e ± 0.325071ei, v* = ±0.650142ei and A* = 
— 0.291667e, clearly indicating the unphysical nature of 
the fixed point. The fixed points are real only when 
to > 25.1. Let us look at the behavior of the flow near the 
disordered p-wave fixed point when we are in this regime. 
For instance, when to = 50, out of the pair of possi- 
ble fixed points there is only one that is physically rele- 
vant. This fixed point is given by u* = 0.0109089e, v* — 
-0.0287184e, A* = 0.136896e. The eigenvalues obtained 
by linearizing around this fixed point for the case to — 50 
are given by -2.49734e, 0.798418e ± 1.57903ei and -e, 
clearly indicating the unstable nature of the fixed point. 
This state of affairs wherein one gets physically relevant 
values for the fixed points that are however unstable con- 
tinues until we reach to ~ 211. For to values that are 
greater than ss 211, we can show that the fixed point 
value of A is negative, thus implying that disordered p- 
wave fixed point goes once again unphysical. 

In conclusion, in this brief report we have studied the 
effect of weak impurities on the first order transition 
found in a p-wave superconductor. It was seen that for 
the physically relevant situation of small number of order- 
parameter components (with to ^ 1), the phase tran- 
sition is still of the fluctuation driven first order type. 
However, it was also shown that the case to = 1, maps 
itself onto the disordered s-wave superconductor that was 
studied in Ref. 0]. 

The conclusion that the critical behavior in a random- 
mass disordered p-wave superconductor remains weakly 
first order is encapsulated in Fig. which depicts the nu- 
merical solutions of the flow equations, Eq. [2j These are 



plotted in Fig. [2 wherein the case to = 1, is contrasted 
with the case to = 2. For the case to = 2, we see the run- 
away flow persists, indicating that the transition remains 
first order. Fig. O clearly shows that for to = L we flow 
into the disordered s-wave fixed point of Ref. 2] . When 
the number of order parameter components is large i.e., 
to > 211, the pure fixed point found in Q goes stable. It 
should be added that our results were obtained in a 4 — e 
expansion. However, in the case of the ordinary s-wave 
superconductors with large values of the Ginzburg pa- 
rameter, k, it has been shown that the 4 — e expansion is 
unable to capture the continuous second order nature of 
the transition seen by using other techniques, [HI], [l6j |. 
Hence, one had to lay recourse to analytical techniques 
other than the 4 — e expansion to capture the continuous 
transition found in superconductors with large k. For in- 
stance, in 1J], the RG is performed in exactly d = 3, to 
capture the continuous phase transition. Efforts are on 
to perform a similar analysis for the case of the p-wave 
superconductor [13]. 
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